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-continued 

trjj{sm) = < Constraint 10.2 

4=1 <=1 



. . . , m(sm), ...,/? 

G(0) represents the nominal model, G + and G"represent the 
allowed upper and lower bound of the model. a y "and a^are 
allowed upper and lower bound for the singular values. 
m(sm) is the dimension of the sub-matrix sm. p is the total 
number of nearly collinear sub-matrices G nc (sm). 

[0056] Equation 10 is the objective function. The objec- 
tive function minimizes deviation from the nominal model 
as long as the required perfect collinearlity can be achieved. 
Constraint 10.1 represents the allowed variation for each 
model element. Constraint 10.2 is contributed from each 
nearly collinear sub-matrix. All u u (sm) and v H (sm) are the 
singular vectors calculated from the original sub-matrix 
G nc (sm), and hence the same directionality is always main- 
tained. 

[0057] In some embodiments, o^.'and a^are set to 

[0058] 1) a i r(sm)=a ij + (sm)=0 if i»j , which corresponds 
to the off-diagonal portion of the singular value matrix; 

[0059] 2) a ii "(sm)«a ii + (sm)=»0 if i>r{sm), which corre- 
sponds to those small singular values to be zeroed out; 

[0060] 3) a u -(sm)o CTi (0)*(l-eps) and a a + (sm)» 
ai(0)*(l+eps) if i <-r(sm), where 0<eps<l is a con- 
stant. Choosing a large value for eps allows large 
variations in the singular values. Since the objective 
function always tries to find the minimal variation for 
the model matrix, it is expected that variation of the 
singular value will also be very mild. Hence, a small 
eps (for instance, eps=0.1) can be safely used. 

[0061] Finally, Equation 10 is a standard QP formula and 
can be solved globally and efficiently. 

Uncollinearization 

[0062] If the process is not collinear, or is nearly collinear 
but needs the controller to fully explore its capability, then 
adjustments to the model can be made to improve the 
condition number so that an improved robustness can be 
achieved. 

[0063] Uncollinearization should satisfy the following 
requirements: 

[0064] 1. The repaired model should have the same 
directionality as the original model because changing 
the direction can cause unwanted control problems that 
can potentially result in a less desirable performance 
than with the original collinearity. 

[0065] 2. The directionality change should be made 
within allowed ranges. Additional restrictions can also 
be imposed, such as additional linear equality or 
inequality constraints. 

[0066] 3. When treating a model matrix larger than 2x2, 
a collinear sub-matrix can share common elements with 
another collinear sub-matrix. Such a case can result in 
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a "zigzag game" or never ending loop, with repairs to 
one sub-matrix giving rise to a need for repairs to the 
other sub-matrix. Hence, the methodology should be 
able to deal with multiple sub-matrices in a synchro- 
nized way. 

To achieve these goals, the following optimization for- 
mula was created: 



Maximize (^^H'/^/t*")^ 0 )) 2 
subject to: 

m(sm) 

gt {sm) = <Yj<r k {sm) * u it {sm) * {sm) <=g~ (sm), 

4-1 

i, J = 1 m(ms) 

07 {sm) = < a- t {sm) <= &• {sm), i = 1 m{sm) 



Equation 11 



Constraint 11.1 



Constraint 1 1 .2 
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£V t (<?) * uik(q) * vjk (<?) = J] <r i (') * urn (0 * vy* (0, 



Constraint 11.3 



for those i j which point to the same element in G, sm, q, 
W , . . . , p, q - t . 



[0067]*'^.^^ 



g'j "denote the upper and lower bounds of 
the allowed model adjustment, cfand a^are upper and 
lower bounds on the singular values, o r+i (0) represents the 
original singular value, and m(sm) is the dimension of the 
sub-matrix sm. p is the total number of nearly collinear 
sub-matrices G nc (sm). Additional explanation of Equation 
11 is provided below. 

[0068] The objective function is to maximize the portion 
of the smaller singular values of all nearly collinear sub- 
matrices. The weighting factor, 1/^^(0), means the 
smaller the original singular value, the more improvement 
the optimizer will attempt to obtain. Constraint 11.1 repre- 
sents the allowed variation for each model element. All 
u^sm) and v jk (sm) (as well as u^q) , v jk (q) , u^t), and 
v jk(0) are tne singular vectors calculated from the original 
sub-matrix G nc (sm) (or G nc (q), G nc (t)), and hence the same 
directionality is always maintained. Constraint 11.2 is the 
allowed singular value variation range, which is discussed 
further below. 

[0069] Constraint 11.3 is needed if there are two sub- 
matrices who share the same element in the original model 
matrix. In this case, any adjustment made for one sub-matrix 
will automatically be coordinated with another associated 
sub-matrix, and hence remove the "zigzag problem." 

[0070] The final goal is to maximize the condition number. 
Equation 11 does not explicitly employ a condition number 
because directly optimizing a condition number will pose a 
computation problem so difficult as to be unrealistic. As 
such, an approximation is made by maximizing the portion 
of the smaller singular values, while keeping the other 
singular values from dropping too low. For this purpose, the 
bounds for each singular value should be set as follows: 



